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Despite 180 years of theory on magnetism, it appears 
that the practice of calculating forces on magnetic media is 
ambiguous, as illustrated by a recent article in this Journal 
[1]. Potentially troubling issues include: Which field B or H 
should be used? Should the total field be used, or only the 
external field? And if the latter, what is meant by "external" ? 
Can/should effects of magnetization currents and/or fictitious 
magnetic poles be included? What is the force on only a por- 
tion of a circuit? We review several well-known approaches 
to magnetic force calculations on elements of rigid circuits, 
and find it very helpful to use an explicit example to compare 
and contrast the various methods. Our discussion reinforces 
that of the most authoritative texts [2,3], but corrects in an 
important way a previous attempt at a systematic review in 
this Journal [4-7] . 

I. INTRODUCTION 

The phenomenon of magnetism was first manifested 
via interactions of bulk magnetic materials (magnets). 
Following the discovery by Oersted [8] that pairs of cur- 
rents exerted forces on one another, Biot and Savart [9] 
identified a corresponding force law, and Ampere [10] 
made the conjecture that all magnetic phenomena are 
actually due to currents, some of which may be bound 
in "molecules" of magnetic materials. This view con- 
trasted with the elegant work of Poisson [11] in which 
forces on magnetic materials (without conduction cur- 
rents) could be explained in terms of a density of mag- 
netic poles bound in the media in such a way that single 
poles can never be isolated. 

Pedagogic discussions of forces on magnetic media tend 
to treat the case of current-carrying conductors sepa- 
rately from that of materials (typically conductors) that 
have nonzero bulk magnetization. The present article 
addresses the general case of magnetic media in which 
conduction currents are flowing, but is restricted to the 
case of forces on elements of a rigid circuit. The litera- 
ture on this topic appears to be both sparse and erratic, 
and we have the (possibly unrealistic) hope of clarifying 
the record. 

Because the details prove to be intricate, the best im- 
pression that a reader may form is simply that one must 
approach magnetic force problems with caution. 

Some more specific conclusions from this effort are: 

1. The consequence of Newton's 1st law that an ob- 
ject does not exert a net force on itself is a guiding 
principle in the analysis of static forces on rigid cir- 
cuits. See sec. III. We do not pursue here the more 
complex case of deformable magnetic media. 

2. In calculating the force on a circuit (or on a portion 
of a rigid circuit) one can use only that magnetic 



field which existed prior to the introduction of the 
circuit, since the fields produced directly (or indi- 
rectly via induced magnetization) cannot cause a 
net force on the circuit (or portion of a rigid cir- 
cuit). See sec. III. 

3. Methods also exist in which the total magnetic field 
can be used in the force calculation, even if the 
force on only a portion of a rigid circuit is de- 
sired. The most straightforward of these uses the 
Maxwell stress tensor (sec. VIII), which the author 
recommends be given more prominence in peda- 
gogic treatments of electrostatic and magnetostatic 
forces. 

4. Considerable care must be given in the treatment 
of systems that include an interface across which 
the magnetic permeability changes discontinuously. 
See sees. IV-IX. 

5. If the Helmholtz bulk force density (calculated from 
the total magnetic fields) is used to find the force by 
integration over a volume that contains such an in- 
terface, the volume integral must be supplemented 
by a surface integral over the (interior) interface. 
See sec. IX. 

6. In a calculation of the force that uses the concept 
of magnetization currents, the magnetic field H 
should be used rather than the field B, and if the 
force is desired on only a circuit element, the field 
must be as existed prior to the introduction of 
the circuit (and not merely the introduction of the 
circuit element). See sec. VI. 

7. In a calculation of the force that uses the concept 
of magnetic poles, the magnetic field B should be 
used rather than the field H, and the prior field B.^ 
should be used if the force on only a circuit element 
is desired. See sec. VII. 

While conclusions 5 and 6 may be known to those 
expert with magnetic media, the author has not found 
them in the literature. Because several of the conclu- 
sions amend conventional pedagogic wisdom, it has been 
critical to verify these conclusions with an example in 
which both I x B forces and magnetization forces are 
present, and which contains an interface with different 
permeabilities on either side. 

II. AN EXAMPLE 

What is the force per unit length on a wire of radius 
a and (relative) permeability fi' when it carries uniform 
conduction current density 



1 



'cond 



7ra^ 



(1) 



and is placed along the z axis in a magnetic field whose 
form is 



Bi = Box + 



X ^ y ^ 
-X- -y 

a a . 



(2) 



before the wire is placed in that field? The medium sur- 
rounding the wire is a liquid with relative permeability 
/X ^ 1. The wire is assumed to be part of a rigid circuit 
that is completely by a loop at "infinity" . 

Note that in asking for the force on the wire along the 
z axis, wo sock the force on only a portion of the complete 
circuit. 

The form of the initial magnetic field has been chosen 
so that there will be both a J x B force associated with 
the uniform field Bq = -Bqx, as well as a force due to 
the interaction of the induced magnetization with the 
nonuniform field Bi = i?i[(.T/a)x — (t//a)y]. 

We have confirmed by the four methods of calculation 
discussed in sees. VI-IX that the force per unit length on 
the wire is, in Gaussian units. 



where the magnetic fields H and 
B = H + 47rM, 



(3) 



(4) 



are related in a linear, isotropic medium of relative per- 
meability \x by 



B = /iH, 

and hence the magnetization density M obeys 



M 



1 



47r 



H. 



(5) 



(6) 



The magnetic field H for this example is given in the 
Appendix, and details of the calculations can be found 
in [12]. 

If the initial magnetic field is uniform {B\ = 0), then 

there is no net force on the magnetization of the wire or 
surrounding medium, and the force per unit length on 
the wire has the very simple form 



IBn 



(uniform external field). 



(7) 



The expression (7) has been verified in a recent ex- 
periment Casperson [1] . See also a related discussion by 
Lowes [13]. 

Strictly speaking, eq. (7) describes the force on the 
conduction electrons, and not that on the lattice of posi- 
tive ions through which the electrons flow. The force (7) 
results in a slight rearrangement of the distribution of the 
conduction electrons and positive lattice ions so that a 
transverse electric field is generated that acts on the lat- 
tice to provide the force experienced by an observer who 
holds the wire at rest. See [14] for further discussion. 



III. THE BIOT-SAVART FORCE LAW 

The result (7) is to be expected from a simple argu- 
ment, consistent with Newton's first law that an object 
does not exert a net force on itself. Namely, the Biot- 
Savart force law (for media of unit permeability) on a 
current-carrying circuit a due to a current-carrying cir- 
cuit b is 

dlaX f 2 = — f dlaX Bj,{a), 



Fa 



where 



'ab 



Bb{a) = 



dh X rob 



(8) 



(9) 



lb 'ab 

In the example of sec. II, the initial field Bj plays the 

role of the field Bf, not due the current in wire a, which 
quickly leads to the result (7). The reader may object 
that the example does not involve media of unit perme- 
ability, so the Biot-Savart iorcc law may have to be mod- 
ified in such a way as to lead to a different result than 
(7). This issue will be pursued in the following section. 

The statement that an object cannot exert a net force 
on itself tacitly presumes that the object has mechani- 
cal integrity as a whole, and will generate internal elas- 
tic forces to counteract possible electromagnetic forces of 
one part of the object on another. If a magnetic circuit 
is mechanically ficxible, we will not be content with an 
analysis of the force on the circuit as a whole. In ad- 
dition, we desire a calculation of magnetic force on an 
element of the circuit, imagining it to be mechanically 
(but not electrically) disconnected from the rest of the 
circuit. In this case, the magnetic field to be used in 
differential form to the Biot-Savart force law is the field 
Bext due to all sources outside the element itself, 

dF = ^dl X Bext. (10) 

The differential expression (10) can have meaning inde- 
pendent of the integral form (8) only if the circuit de- 
forms, in which case the problem is not one of statics. 
It is well known that the differential force (10) does not 
satisfy Newton's 3rd law in cases of isolated current el- 
ements, because in such dynamical systems the electro- 
magnetic field momentum is varying [16]. Newton's 3rd 
law is respected via the appropriate electromagnetic ver- 
sion of his 2nd law. 



E 



dP 



mech 



dP 



EM 



dt 



dt 



= 0, 



(11) 



where Pmech and Pem are the mechanical and electro- 
magnetic momenta, respectively. For a discussion of dc- 
formable circuits, see [17]. The topic of forces on mag- 
netic liquids, including striction effects, is treated in [18]. 

In the remainder of this article we assume that the 
object on which we desire to calculate the magnetic force 
is a rigid body at rest. 
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IV. MICROSCOPIC AND MACROSCOPIC 
MAGNETIC FIELDS 

In the macroscopic approach to calculation of forces on 
magnetic media one considers in general a current density 
J (electric charge crossing a directed unit area per unit 
time) and a net magnetization density M (net magnetic 
dipole moment per unit volume) as well as the magnetic 
fields B and H = B — AttM. In a microscopic view one 
considers collections of individual charges and/or mag- 
netic moments, but averaged quantities like the magneti- 
zation density M arc not yet defined. Hence, the fields B 
and H are identical in the microworld, and it is a matter 
of convention which symbol is used for the microscopic 
magnetic field. 

The symbol H was used for the magnetic field by the 
early workers Ampere, Biot and Savart, and Poisson [11], 
which led authors such as Lorentz [19] and Landau [3] to 
use this symbol for the microscopic magnetic field. The 
symbol B was introduced by Thomson around 1850 [21] 
in the form (5), which suggests that it is to be derived 
from the more fundamental (or anyway more familiar at 
the time) field H. However, as apparently first noted by 
Lorentz [19], the macroscopic average of the microscopic 
magnetic field is B and not H. 

To maximize the continuity between the microscopic 
and macroscopic views, the author prefers that the mi- 
croscopic magnetic field be labeled B. Then, the force 
on an electric charge e with velocity v in microscopic 
electromagnetic fields E and B is 



(12) 



the Lorentz force law, which also represents the average 
force on the charge in the presence of macroscopically 
averaged fields E and B. It is always understood that the 
fields E and B in the Lorentz force law do not include 
the fields of the moving charge itself. 

High-energy physicists such as the present author, con- 
sider the Lorentz force law (12) to be experimentally con- 
firmed (and continually reconfirmed) for over 50 years 
[22] in the case of high-energy particles moving inside 
media such as iron where B = /iH ^ H. That is, when 
applying the Lorentz force law a charge inside a macro- 
scopic medium, the appropriate macroscopic average of 
the microscopic magnetic field is indeed the macroscopic 
field B. This insight is affirmed in sec. 8.2 of [15], and 
sec. 22.1.1 of [23], 

We therefore expect that other methods of calculating 
forces on currents in macroscopic media will be consistent 
with the Lorentz force law using the macroscopic field B. 
In particular, we expect that the Biot-Savart force law 
for a current I in a macroscopic magnetic field Bj not 
due to that current itself is 



F = - I Idl X Bi = - (f) dl X Bi, 



(13) 



no matter what permeabilities exist. Thus, we reaffirm 
that eq. (7) is the force on the wire in our example, if the 
initial field is uniform. 



V. THE BIOT-SAVART FORCE LAW IN A 
PERMEABLE MEDIUM 

Despite the simplicity of the result (13), care is needed 
when using the Biot-Savart force law in permeable me- 
dia. We review this issue by starting with the case that 
all wires and their surrounding media have the same per- 
meability jjL ^ 1. Then there is neither a surface current 
nor a fictitious pole density at the interface between the 
wire and the liquid. However, there remains a volume 
current density 



cV X M 



1 



A-k/c 



V X H = (m - l)Jcond, (14) 



using Ampere's law. 



47r 

V X H = Jcondj 

C 



(15) 



that relates the magnetic field H to the conduction cur- 
rent density Jcond- Thus, the total current density is 



Jtotal — M Jcond- 



(16) 



The fact that the total current density (16) does not 
equal the conduction current density in a permeable 
medium contradicts the view of Lorentz [19], as reaf- 
firmed in sec. 30 of [3]. This appears to be one of the 
very few oversights in these distinguished works. 

The force (13) on the wire that carries conduction cur- 
rent density Jcond can be written as 



Jcond X Bj cfVol. 



(17) 



If instead we wish to use the total current density (16) 
we must write 



F = - / xBidMo\=- / Jtotal X H, dVol. (18) 

H c ' 



Another aspect of the analysis of Biot and Savart is the 
calculation of the magnetic field from the current density. 
The microscopic version of Ampere's law. 



An 

V X B = Jtotal, 

c 

corresponds to the prescription that 

Jcond X ^ 



(19) 



g ^ 1 ^ JtotalXf 



-J 



cfVbl = yuH, 

(20) 



supposing the permeability is uniform and may be taken 
outside the integral. Hence, the macroscopic version of 
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Ampere's law, eq. (15), corresponds to the prescription 

that 

H = i / J^o-^d ^^ rfvol = V X i / ^dVol, (21) 
c J r'^ c J r 

independent of the permeabiHty. This result, combined 

with Ampere's law (15), is consistent with Hclmholtz' 
theorem [15] provided V • H = 0, as holds within a 
medium of uniform permeability. 

The form of the cq. (8) for the force on circuit a due to 
circuit b supposing the wires and the surrounding media 
all have permeability n is therefore 

Fa = f dlaX f 2 ' (22) 

Ja Jb '''ab 

where and /(, are the conduction currents in the cir- 
cuits. 



and the total current density is Jtotai = /i'/z/7ra^ from 
cqs. (1) and (16). Then, evaluation of eq. (24) leads to 
cq. (3). 

The correct result for the force on a portion of a rigid 
circuit would not be obtained from eq. (24) if we used the 
initial magnetic field Bj, or if we used the total magnetic 
fields B or H on the wire. Furthermore, we would not 
obtain a correct result if we used the field in a vacuum 
cavity of radius a at the position of the wire cither before 
or after the wire was inserted. The proper initial field 
Hj is the one before the wire was inserted into the liquid 
dielectric. 

When ^ /k' ^ 1, we cannot rewrite the first term of 
eq. (24) as / (Jcond/c) x Bj cTVol, which would incorrectly 
suggest that we could ignore the volume magnetization 
current density Jm but not the surface current density 
Km. 



VI. THE BIOT-SAVART FORCE LAW PLUS 
BOUND CURRENT DENSITIES 

We also see that eq. (22) holds even if the wires have 
permeabilities /Xo and that diflfer from the permeability 
/i of the surrounding medium, since the magnetic field 
due to wire h at the position of wire a before wire a 
was introduced is given by Bt = /iHb, which depends on 
neither /ia nor /z;,. However, in this case there will exist 
effective surface currents, 



K 



M 



cAM X n, 



(23) 



at the interface between the wires and the surrounding 
medium, where AM is the difference between the mag- 
netization on the two sides of the interface, and n is the 
unit normal to the interface. The force on these surface 
currents follows the form (18) that uses the magnetic field 
H rather than B. Hence the total force on a permeable 
wire surrounded by a permeable medium can be written 
as 



c 



^ j Jtotai X Hi cTVol + i y" Km X Hj dS. 



(24) 



The is the appropriate version of the Biot-Savart law 
if we wish to include magnetization forces via the so- 
called bound current densities. However, the Coulomb 
Committee in their eq. (1.3-4') [6], and Jefimenko in his 
eq. (14-9. 13a, b) [20], recommends that the initial field 
Bj be used rather than Hj, which would imply a force /U 
times the above. 

For the example of sec. II, the surface current density 
is obtained from eq. (6) and the total magnetic field H, 
eq. (43), as 



Km = 



47r 



21 



2/xc 



II + 11' 



-{Hosm.6 + H-isva.2e) 



(25) 



VII. THE BIOT-SAVART FORCE LAW PLUS 
FICTITIOUS MAGNETIC POLES 

Following Poisson [11], the forces on the magnetization 
of the media can also considered as due to a density of 
fictitious magnetic poles, rather than being due to cur- 
rents Jm and Km- Some care is required to use this 
approach, since a true magnetic pole density pM would 
imply V • B = 47rpM, and the bulk force density on 
these poles would be F = pmB. However, in reality 
= V- B = V-(H-h AttM), so we write 



V • H = -47rV • M = 47rp 



M, 



(26) 



and we identify pu = — V • M as the volume density of 
fictitious magnetic poles. Inside linear magnetic media, 
such as those considered here, B = /iH and V ■ B = 
together imply that pM = 0. However, a surface density 
(Tm of fictitious poles can exist on an interface between 
two media, and we see that Gauss' law for the field H 
implies that 



(H2 - Hi) ■ n 

47r 



(27) 



where unit normal fi points across the interface from 

medium 1 to medium 2. The surface pole density can 
also be written in terms of the magnetization M = 
(B - H)/47r as 



(Tm = (Ml - M2) • n. 



(28) 



since V • B = insures that the normal component of B 

is continuous at the interface. 

The force on the surface density of fictitious magnetic 
poles is 



(TMBj 



(29) 



since the fictitious poles couple to the macroscopic aver- 
age of the microscopic magnetic field, as anticipated by 
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Thomson and Maxwell [24]. Equation (29) is in agree- 
ment with prob. 5.20 of [26]. However, the Coulomb 
Committee in their eq. (1.3-4) [6], and Jefimenko in his 
eq. (14-9. 9a, b) [20], recommends that the initial field Hj 
be used rather than Bj when using the method of fic- 
titious magnetic poles, which would imply a force l/fi 
times the above. 

The total force on the medium in this view is the sum 
of the force on the conduction current plus the force on 
the fictitious surface poles, where to avoid calculating a 
spurious force of the rigid wire on itself we use the initial 
magnetic field B,, 



Jcond X Bj cfVol -I- / CTMBj dS. 



(30) 



In the example of sec. II, the density of fictitious mag- 
netic poles on the surface r = a is given by 



Cm = 



Hr{r = a+) - Hr{r = a") 



An 



1 n-n' 



2-n jj, -\- ji' 



- {Ho cos 6 + Hi cos 29). 



(31) 



Then, evaluation of eq. (30) again leads to the result 
(3). This reaffirms that one should use the initial field 
Bi and not Hi and not the total field B in the Biot- 
Savart method (30) with fictitious magnetic poles when 
calculating the force on only a portion of a rigid circuit. 



VIII. THE MAXWELL STRESS TENSOR 

The methods of calculating the force on a circuit dis- 
cussed thus far are require care in that only the fields Bj 

and Hj prior to the addition of the circuit appear explic- 
itly in the force calculation, although the total fields B 
and H are needed in the calculation of the magnetization. 
The erratic literature on this topic is ample evidence that 
confusion as how to implement these calculation is likely. 

Hence, it may be preferable to use methods that in- 
volve only the total fields B and H. In the author's view, 
the most reliable general method for the calculation of 
electromagnetic forces is that based on the Maxwell stress 
tensor [29] , which is a formal transcription of Faraday's 
"tubes of force". The j component of the electromag- 
netic force F on the interior of a closed volume in a linear 
medium with oriented surface element dS is given by 



dSk, 



(32) 



where 



1 S 
Tjk = ^EjDk + BjHk - -^(E • D + B ■ H). (33) 



The form (33) ignores interesting strictive effects in com- 
pressible media [28]. 



The surface over which the stress tensor is integrated 
need not correspond to a physical surface, which leads to 
the question of how the electromagnetic force is trans- 
mitted to the physical matter inside that surface. The 
answer is, of course, via the electromagnetic fields that 
enter into the stress tensor. Prior to, and even some- 
what after. Maxwell, more physical explanations were 
considered necessary, which led to the diversionary search 
for the aether that, among other activities, would trans- 
mit the forces from the imaginary surface to the matter 
within. 

We have already noted in sec. II that even the J x B 
force concept is at least once removed from a force on 
the positive ion lattice of a conductor - which caused 
confusion to Maxwell [30]. 

For the example of sec. II, we calculate the force on 
unit length of the wire by integrating the Maxwell stress 
tensor over a cylindrical surface of radius r > a, so that 
any effects at the surface r = a are included. The result 
is, as expected, given by eq. (3), independent of radius r. 

Equation (3) for the case that Bi = was deduced in 
a similar manner in ref. [2]. 

If we integrate the stress tensor over a cylinder of ra- 
dius r < a the result is 



c 



(34) 



Since the limit of this as r ^ a does not equal the result 

(3) for r > a, we infer that there are important effects at 
the interface r = a. The permeable liquid is presumably 
contained in a tank of some characteristic radial scale 
6 ^ a, at whose surface additional magnetization forces 
will arise. We consider these forces as distinct from those 
at the interface r — a, and that only the latter are part 
of the forces on the wire. 



IX. THE HELMHOLTZ BULK FORCE DENSITY 

An expression for a bulk force density f in magnetic 
media can be obtained by transformation of the surface 
integral of the stress tensor into a volume integral. See, 
for example, sees. 15 and 35 of [3]. The result, again 
ignoring magnetostriction, is 



1 H 
f = -Jcond X B - — Vm, 

c Stt 



(35) 



which is due to Helmholtz [27] . As for the Maxwell stress 
tensor, the fields B and H in eq. (35) are the total fields 
from all sources. 

However, eq. (35) is not sufficient for the case that 

the permeability takes a discontinuous step at an inter- 
face within the volume of interest. To see this, we recall 
the usual derivation of the bulk force density, beginning 
with Gauss' law to transform the surface integral of the 
Maxwell stress tensor into a volume integral of a force 
density f. 
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Fi 



Tij dSj 



fi dVol. 



(36) 



Recalling eq. (33), the force density f is given by 



f = VT = 



1 

47r 



(B • V)H + H(V • B) - iv(B • H 



1 



(37) 



since V-B = always, and we first suppose that B = /xH 
involves a continuously varying permeability. The usual 
argument then proceeds by noting that 



Vil^ = 2(H • V)H + 2H X (V X H) 

= 2(H • V)H - — J,o„d X H, 
c 



(38) 



using Ampere's law (15). Inserting cq. (38) in (37) we 
arrive at eq. (35). However, if the volume of interest in- 
cludes an interior interface across which the permeability 
takes a discontinuous step, wc should revert to the first 
form of eq. (37) when performing the volume integral 
across the interface. Defining n to be the unit normal to 
the interface, and noting that B • H = ^Hf + B^/ij, where 
the tangential and normal components, Ht and are 
continuous across such an interface, the resulting surface 
integral is 



1 

4tt 



BnAHn - -A(B ■ H) 



ndS, 



The total magnetic force acan now be written 



(39) 



+ 



cond X B dVol / iJ^V/z dVol 

c 8tt J 

l-i.{^)fBln4S-^jH^u,S (40) 



where the surface integral is over the interface which is, 
in general, interior to the volume of integration, and AA 
is the difference of quantity A on the two sides of the in- 
terface. To apply eq. (40) to example of sec. II, we note 
that the second integral is zero (considering the integral 
over the interior interface as distinct from the bulk vol- 
ume integrals), and that n = f on the interface at r = a. 
Therefore, wc need only those terms of Bf. and Hg that 
vary as sm.0 or cos 9. Referring to the Appendix, we find 



A ( i j Bl - AfxH^ 

^IJ,- H' ( IBq 



■ sm( 



Bo Hi 



■ cos I 



+ fj,' \ ca 2 
Using this in eq. (40), we again obtain the result (3). 



+ ... (41) 



X. APPENDIX: THE FIELD H FOR THE 
EXAMPLE 

The magnetic field H for the example of sec. II is, in 
both cylindrical and rectangular coordinates. 



H= i 



■j^{H,cose + H,r.cos2e)i 



1 



At—A*' o.^ 



(r < a). 



cos f 



Hi 



f^ — fJ' a' 



cos 20 



-Hi {^ + ^^)siii2e\0 

-'^^^(^+^{H, + Hilcose)\± 
+ {^^ose-j^Hi^-sme)y 
sine + Ho fl- ^4 cos 26(1 



(r > a), 
(r < a), 



+ 
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- Hi f r cos 6' - 4 cos 3611 
cos6'- ^i^ifo4sin26' 



Hi 



(^^sinO 



+ ^^sin3e)]y (r>a). 
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